Polar Charge Fluctuation and Superconductivity in Organic Conductor 
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Superconductivity and polar charge fluctuation are studied in an organic conductor with the dimer-molecule 
degree of freedom. The extended Hubbard models, where the intra-dimer electronic structure and the inter- 
dimer Coulomb interactions are taken into account, are analyzed by the random-phase approximation and the 
fluctuation-exchange approximation. Superconductivity appears in a vicinity of the charge-density wave (CDW) 
phase where the electronic charge distributions are polarized inside of dimers. The extended s-wave type paring 
is favored and its competitive relations with the superconductivity due to the spin fluctuation depends on the 
triangular lattice geometry. Comparison between two superconductivities realized near the polar and non-polar 
CDW phases are also presented. 

PACS numbers: 74.70.Kn, 77.80.-e, 74.20.Mn 



Unconventional superconductivity is one of the central 
issues in modern solid state physics. Beyond the stan- 
dard Bardeen-Cooper-Schrieffer theory, non-phononic mech- 
anisms and anisotropic pairing symmetries have been pro- 
posed theoretically and examined experimentally in a wide va- 
riety of materials, such as copper oxides, iron pnictides, heavy 
fermion compounds and so on. Organic conductors are one 
of the families in which exotic superconductivities have been 
examined intensively and extensively since the discovery of 
superconductivity in the TMTSF salts^ 

Recently much attention has been attracted in layered 
quasi two-dimensional K-type BEDT-TTF salts, k-(BEDT- 
TTF)2X (X: a monovalent anion)^ The maximum super- 
conducting (SC) transition temperature (T c ) of 13. 2K is ob- 
served in k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 under high pres- 
sure. 5 In the conducting BEDT-TTF layers, two BEDT-TTF 
molecules are regarded as a molecular dimer which builds a 
two-dimensional triangular lattice. A geometrical anisotropy 
in the triangular lattice is adjustable by changing an anion 
X and/or applying pressure. One hole-carrier per dimer oc- 
cupies the lowest-unoccupied molecular-orbital in a dimer. 
When the intra-dimer Coulomb interaction is much stronger 
than the band width, the system is regarded as a Mott insu- 
lator, which is termed a so-called dimer-Mott insulator^ 7 . An 
antiferromagnetic ordered state as well as a spin-liquid state 
observed in a series of materials are consequences of the Mott 
insulating state£^ Since the SC phase appears at a vicinity 
of the magnetic phases^ a mechanism of the superconductiv- 
ity has been often examined from the view point of the spin 
fluctuation^^ 

Recently, dielectric anomaly is reported in one of the K-type 
salts, k-(BEDT-TTF) 2 Cu 2 (CN) 3 ^ A cusp-like structure and 
strong frequency dependence are observed in the temperature 
dependence of the dielectric constant around 30K. One possi- 
ble origin of this dielectric anomaly is attributed to the electric 
dipole moments inside of the BEDT-TTF molecular dimers, 
i.e. an electronic-charge distribution without the space- 
inversion symmetry in a dimer. 18 ~ 2 - ) Similar dielectric anoma- 
lies are also reported in k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C^ and 
/3'-(BEDT-TTF) 2 ICl 2 ifii These experimental results trigger 
reinvestigations of the electronic structure in the dimer-Mott 



insulating systems, and a mechanism of the superconductivity. 

In this paper, we study superconductivity in low- 
dimensional organic conductors, in particular, focus on a role 
of the polar charge fluctuation. We introduce the extended 
Hubbard model where, in contrast to the previous works, we 
take into account both the electronic degree of freedom in- 
side of a dimer and the inter-dimer Coulomb interactions, 
which are essential for the polar charge order and fluctua- 
tion. Two types of the molecule configurations where two 
and four molecules are included in a unit cell are analyzed by 
the random-phase approximation (RPA) and the fluctuation- 
exchange (FLEX) methods. The SC phase appears near 
the charge-density wave (CDW) phase where the electronic 
charge distribution is polarized inside of dimers, as well as 
near the spin-density wave (SDW) phase. The extended s- 
wave symmetry pairing is favored near this CDW phase and 
is realized cooperatively with the d xy -type pairing due to the 
spin fluctuation. Contribution of the polar charge fluctuation 
is remarkable around the so-called frustration point where the 
dimer molecules are arranged on an equilateral triangular 1 at - 
tice. A reentrant feature of the phase boundary is observed for 
the superconductivity induced by the polar charge fluctuation. 

We introduce the two types of the molecule configurations, 
termed the two-band model and the four-band model, where 
inequivalent two molecules and four molecules are introduced 
in a unit cell, respectively. Schematic views are shown in 
Fig. Q] First we present the Hamiltonian for the two-band 
model: 
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where Ci aa is an electron annihilation operator at the i-th 
dimer with molecular orbital j(= a,b), and spin cr(=t, ■!•)> 
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tor. The first term represents the intra-dimer electron tran sfer, 
and the second and third terms represent the intra-molecule 
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the bubble-type diagrams are only considered, for simplic- 
ity The spin susceptibility x s q {iuj n ) and the charge sus- 
ceptibility Xq(iw n ) m the RPA scheme are given by 



(3) 



where — and + signs are for the spin and charge susceptibili- 
ties, respectively, "hat" implies the matrix, and x q {iu n ) is the 
bare susceptibility with the Matsubara frequency u> n defined 
by 



FIG. 1: (Color online) Schematic views of the two-band model in 
(a) and the four-band model in (b). The Coulomb interactions, V, 
and V' in the two band model, and V&, V v 



model are introduced in the same ways as t and t', and tb, t p 
respectively. Ovals represent the BEDT-TTF molecules. 



and t„ 



and inter-molecule electron-electron interactions inside of a 
dimer, respectively. The last two terms describe the inter- 
dimer electron transfers and electron-electron interactions, re- 
spectively. In the two-band model, the inter-dimer trans- 
fers and Coulomb interactions are assumed to be diagonal as 

£[■ = Sryry/tij and V^ 1 = S^yVij, for simplicity. The ge- 
ometrical anisotropy in a triangular lattice is represented by 
the two kinds of ty, t, t', and those of Vy, V, V , as shown 
in Fig. [Ha). Spatial difference between the two molecules in 
the same dimer is neglected, for simplicity. The average elec- 
tron number per dimer is fixed to be 3, corresponding to the 
averaged electron density in K-type BEDT-TTF salts. The in- 
teraction terms, i.e. the U, U' and V^ 1 terms in Eq. (Q~|i, are 
represented in a unified fashion as 
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where Ck/Sa is the Fourier transform of c^, and the interac- 



tions \W° 



are the 4x4 matrices. 



q J(a/3)( 7< 5) 

In the four band model, inequivalent two dimers in a unit 
cell are identified by an index £(= 1,2) and the annihila- 
tion operator introduced in Eq. ([TJ is replaced as c, Qcr — > 
Ci^aa where (£,a) = {(1, a)(l, b)(2, a)(2, &)}. The intra- 
dimer electron transfer, the intra-molecule and inter-molecule 
electron-electron interactions inside of a dimer are the sim- 
ilar forms with the first three terms in Eq. (UJ, respectively. 
Three kinds of the inter-dimer electron transfers, tf,,t p ,t q , 
and the inter-dimer Coulomb interactions, Vb,V p ,V q , are in- 
troduced, as shown in Fig. QIb). The effective interactions 
\W q a ]( 5/ §)( 7 ^ with a = (£, a) introduced in Eq. © are the 
16 x 16 matrices. 

The multi-band Hubbard models with the long-range Coul- 
omb interactions are analyzed by the RPA and FLEX meth- 
odsi22r— We consider the so-called bubble-type and ladder- 
type diagrams where the Coulomb interactions are considered 
up to the infinite order of the perturbation. For the Coulomb 
interactions U and U' , both the ladder- and bubble-type dia- 
grams are adopted, and for the inter-dimer interaction , 
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We define U s q = - W^, \J c q = + W^, and the bare 
Green's function G Q k (ie n ). In the FLEX method, the spin and 
charge susceptibilities are calculated in the equations given in 
Eqs. (0 and where the bare Green's functions are replaced 
by the full Green's functions defined by the Dyson equation, 
G fc (ie„) _1 = G^(ze n ) _1 -£ k (ie n ). The self energy is given 
by 
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where we introduce the effective interaction: 

v q (iu n ) =\u s q x q ^ n )u s q + l -u c q r q {^ n )u c q 
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It is convenient to introduce the following static suscepti- 
bilities: the static spin susceptibility 
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the static charge susceptibility 
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and the static polarization susceptibility 
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° rU Oer TH , 

the Fourier transforms of the spin, charge and polariza- 
tion operators, respectively. These are defined by S" = 

( 1 / 2 )Ea(=a,6)( n »t -niai), Qi = (1/2) Ea(=o,6)<r n ™° 

and P l = 0-/2)J2A n ™<* ~ n «°0 for the two-band 



where 0{r) = e Tfi Oer TH , and S z q , Q q and P q are 
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model, and Sf = (1/2) Ee(=i, 2 )a(=o,6)( n ie«t " »*ifa4.)> 
<5» = ( 1 / 2 )E Q (=a,fc)<r(^iaT - n i2ao -) and Pi = 
(I/ 2 ) S5(=i,2) C r( ri i«ao--^6 C r) for the four-band model. The 
static susceptibilities are represented by the spin- and charge- 
susceptibility matrices previously introduced. The explicit 
representations, for example for the two band model, are given 
as 

x s r = \ E ^ (7171 ) (7272 )(o), do) 

7i/Y2(=a,&) 

charge _ _ ST^ c (n\ , i i % 

Aq 2 A q(7l7l)(7272) V u /' '•'V 

7ii72(=o,6) 

and 

Xq° lar = ^ E ^q(7l7l)(7272)(°) e 7l72; ( 12 ) 
7l:72(=a,6) 

where e 7172 = 1 for -ji = j 2 an d — 1 for 71 7^ 72- The formu- 
lae for the four-band model are given in the similar ways. On 
the analogy of the definitions of the static-susceptibilities, we 
introduce the "bare" static charge-susceptibility x q charge , and 
the "bare" static polarization susceptibility x q polar . These are 
defined in Eqs. (0 and ©, where the expectations (• • • } are 
replaced by those without the Coulomb interactions, and are 
calculated by Eqs. (fTTl i and dT2l . where x q (0) is replaced by 
Xq(0) in the two-band model. 

The SC transition temperature and the gap function are cal- 
culated by the Eliashberg equation. The linearized equation is 
given by 

xG^( AlAl /)(ie I „)G„ fc /(„^)(-ie m )A^ ) ( Ai , l ,, ) (je m ), 

(13) 

where I — s(t) for a spin singlet (triplet) paring. A parameter 
A is introduced and T c is determined by the condition A = 1. 
We define the interactions 

f ¥ = \u s q x s q u s q - \u c q x%u c q + \ (u q + u c q ) , (14) 

and 

^ = ~\u q x q u q - \u q x c q u c q + \ {-u s q + u c q ) ■ (15) 

We take 32x32 (64 x 64) momentum-point meshes in the Bril- 
louin zone and up to 8,192 (16,384) Matsubara frequencies in 
the numerical calculations by RPA (FLEX). 

First we show the numerical results in the two-band model. 
We take t'/t = V' /V > 1 that correspond to the ge- 
ometrical anisotropy in a triangular lattice for k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 . The phase diagram at T/t = 0.01 cal- 
culated by RPA is presented in Fig. EJa). The topologically 
similar phase diagrams are reported in Refs. |29|and|30l where 




kx/n K/tt K/n 



FIG. 2: (Color online) (a) Phase diagram at T/t = 0.01 for the two- 
band model calculated by RPA. (b) A schematic charge configuration 
in the polar CDW phases. Contour maps of the SC gap functions at 
(t/t ,V/t ) = (0.7,0.5), (0.75,0.8) and (0.875,0.86), which are 
marked by circles in (a), are shown in (c), (d) and (e), respectively. 
Bold lines represent the Fermi surface. Parameter values are chosen 
to be U/t = 3, U'/to = 2.1, and t'/t = V'/V = 1.05. 



one molecule without the dimer degree of freedom is intro- 
duced at each site in triangle and square lattices. 

The boundary for the SDW (CDW) phase is identified by 

the condition det[l - (+)x°(0)£>q (c) ] = 0. The SDW and 
CDW phases appear in small t/to and small V/to regions, re- 
spectively. The phase boundary for SDW does not depend 
on V, because the inter-dimer Coulomb interactions are not 
included in The momenta where the static spin sus- 

ceptibility x q pm takes the maxima are q = (0.757T, 0.3757r) 
and (-0.757T, -0.375tt) at (t/t ,V/t ) = (0.65,0.5). To 
examine the CDW phase in more detail, we calculate the 
static charge and polar susceptibilities. At the boundary for 
the CDW phase, x q ° lar diverges and x q harge does not. That 
is, the charge disproportionation occurs inside of dimers and 
the electron number per dimer is equal for all dimers. A 
schematic charge configuration is shown in Fig. [2] The mo- 
menta where Xq° lar takes the maxima are q = (2tt/3, 2tt/3) 
and (-2tt/3, -2tt/3) at (t/t , V/t ) = (0.75, 0.8), implying 
a three-fold periodicity of the electric dipole moment and no 
macroscopic polarization. 

The SC phase boundary and the SC gap functions are shown 
in Fig. |2] The SC phase volume near the SDW phase is larger 
than that near the CDW phase. In the present lattice struc- 
ture with the D211 point group, a spin-singlet SC-gap sym- 
metry is classified by the Ai g and Bi g irreducible represen- 
tations. The dry-type Ai ff SC gap is realized near the SDW 
phase at (t/t ,V/t ) = (0.7, 0.5) and (0.75, 0.8), and is grad- 
ually changed into the extended s-type Ai g gap near the CDW 
phase. This extended s-type gap is also obtained by the cal- 
culation where we set U s x s U s — in Eqs. ( fT4b and dT5b . 
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FIG. 3: (Color online) (a) Superconducting transition temperatures 
near the CDW phase at t'/t = V'/V = 1.05. (b) A ratio of 
the transition temperatures at V/to — 1.5 to those at V = 0, 
i.e. r T = T c (V/t = 1.5)/T C (V = 0) as functions of the 
anisotropy in the electron transfer, t'/t. The transition temperatures 
at V/to = 1-5 are also shown. The charge and spin susceptibili- 
ties at T/t = 0.0025 for t'/t = 0.8 and those for t'/t = 1.05 
are presented in (c) and (d), respectively. All data are obtained in 
the two-band model analyzed by FLEX. Other parameter values are 
chosen to be U /to = 4, U' /to = 2.8, and t/t = 0.5. 



It is concluded that the polar charge fluctuation induces the 
extended s-type SC state. 

The two-band model is further analyzed by the FLEX 
method. The CDW phase with the electric dipole moments 
inside of dimers, obtained by RPA [see Fig. [2), is reproduced, 
but the SDW phase is not realized in finite temperatures due 
to the Mermin- Wagner's theorem satisfied in FLEX^i The 
SC transition temperatures near the CDW phase are shown 
in Fig. |3a). The SC transition temperatures are determined 
by fitting the numerical data of A calculated down to 0.002io 
by a function A = — alogT + b. The transition temperature 
monotonically increases, when the system approaches to the 
phase boundary. This is evidence that the SC state around this 
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FIG. 4: (Color online) (a) Phase diagram in the Vb-V q plane, and 
(b) that in the at-Vt plane at T '/to = 0.1 for the four-band model 
calculated by RPA. Schematic polar and non-polar CDW structures 
are also shown. We choose t^/to = 0.5, t p /to = 0.35, and t q /to = 
—0.13 in (a), and h/to = 0.5a t , t p /to = 0.35a t , and t q /to = 
— 0.13at in (b). Other parameter values are chosen to be U/to = 
1.6, U' /to = 1, and V p /t = 0.95. 



region is caused by the polar charge fluctuation. A numerical 
value of T c for t'/t = 1.05 (t'/t = 0.8) is estimated to be 
about 0.5K (5K) at V/t = 1.5, when< = 0.2eV is assumed. 
It is worth to note that, in comparison with the results by RPA 
[see Fig. |2j a)], the polar CDW phase and the SC phase near 
the CDW phase are realized until at least t/to = 0.3. This 
originates from a suppression of the SDW instability due to 
the Mermn- Wagner's theorem and is consistent with K-type 
BEDT-TTF salts^ 

Effects of the triangular-lattice geometry on the SC tran- 
sition temperature is examined by calculating a ratio of the 
transition temperatures at V/to = 1.5 to that at V = 0, 
i.e. r T = T c (V/t = 1.5)/T C (V = 0) as functions of the 
anisotropy in the electron transfer, t'/t [see Fig. Ob)]. We 
also plot the t'/t dependence of T c at V/to = 1.5, which is 
consistent with the results in Ref. [TH It is clearly shown that 
tt takes its maximum around the so-called frustration point 
t'/t = 1. That is, a relative contribution of the polar charge 
fluctuation to the spin fluctuation increases around t'/t = 1. 
In Figs.Oc) and (d), the spin susceptibility and the polariza- 
tion susceptibility at T/t = 0.0025 are shown for t'/t = 0.8 
and 1.05, respectively. We define K s and K p are the maxima 
of Xq P ' n an d Xq° lar ' respectively, when the momenta q are var- 
ied. It is shown that K s at t'/t = 1.05 is reduced from that in 
t'/t = 0.8, on the other hand, K p are almost independent of 
t'/t. This is because the spin fluctuation is almost governed by 



Xq which is sensitive to the anisotropy of the electron trans- 
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FIG. 6: (Color online) (a) Phase diagram at T/t = 0.01 for the 
two-band model for t'/t = V /V = 0.9 calculated by RPA. Con- 
tour maps of the SC gap function at (t/io, V/to) = (0.6, 0.73), and 
(0.74,0.93), which are marked by circles in (a), are shown in (b) 
and (c), respectively. Bold lines represent the Fermi surface. The 
parameter values are chosen to be U /to = 2.4, and U' /to = 1.6, 



FIG. 5: (Color online) (a) Finite temperature phase diagram for the 
four-band model calculated by RPA. Contour maps of y max (q), the 
largest eigenvalues of the Fourier transform of V^ 1 , are plotted in 
(b) for Vq/to — 0.35 and in (c) for Vq/to = 0.85. Contour maps 
of the "bare" static polar-susceptibility Xq polar and the "bare" static 
charge-susceptibility x% chargc at T/t = 0.05 are plotted in (d) and 
(e), respectively. The parameter values are chosen to be U/to = 1.6, 
U'/to = 1, V p /t = 0.95, Vb/to = 0.9, U/to = 0.5, t p /t = 0.35, 
and tq/to — —0.13. 



fer, t'/t, but the charge fluctuation is dominated by both the 
inter- molecule and inter-dimer Coulomb interactions. 

The SC state induced by the polar charge fluctuation is stud- 
ied in the realistic four-band model. In Fig. Ufa), the phase dia- 
gram in the four-band model calculated by RPA at T/Iq =0.1 
is presented in the Vb-V q plane. Numerical values for the 
transfer integrals are chosen by the extended Huckel tight- 
binding band calculations^ Two types of the CDW phases 
appear; the non-polar CDW for large V q where electron num- 
ber per dimer are different in two inequivalent dimers, and the 
polar CDW for large VJ, where the polar charge distribution 
is realized inside of dimers. Schematic charge configurations 
are shown in Fig. [4] The SC phases are observed near both the 
phase boundaries for the polar and non-polar CDW phases. 
In Fig. SJb), we plot the phase diagram in the at-V> plane, 
where the inter-dimer transfers are given by a parameter at 
as tb/to = 0.5a t , t p /t — 0.35a t , and t q /t = — 0.13a t . 
This phase diagram qualitatively reproduces the results in the 
two-band model (see Fig. [2]). We confirm that the SC phase 



induced by the polar charge fluctuation does not depend on a 
detailed lattice structure. 

Finite temperature phase diagram near the polar and non- 
polar CDW phases calculated by RPA is presented in Fig. [5] 
as a function of V q . A reentrant feature is observed in the 
phase boundary between the polar CDW phase and the metal- 
lic or SC phase, but not in the boundary for the non-polar 
CDW phase. This is interpreted by the competitive momen- 
tum dependences of the charge susceptibility and the effec- 
tive interaction as follows. In the RPA scheme, the momen- 
tum dependences of the susceptibilities are governed by those 
of the bare susceptibilities and the long-range Coulomb in- 
teractions. Contour maps of V -max (qr), which is defined as 
the largest eigenvalues of the Fourier transforms of the inter- 
dimer interaction V^ 1 , are calculated. Results near the polar 
CDW phase (V q /t = 0.35) and the non-polar CDW phase 
{Vq/to = 0.85) are plotted in Figs.[3Jb) and (c), respectively. 
We also show the contour maps of the "bare" static polar- and 
charge-susceptibilities, \q polar an d Xq charg0 , in Fig.|5jd) and 
(e), respectively. It is noticed that momentum dependences of 
X° q polar and y max (q) at V q = 0.35 [see FigsHb) and (d)] are 
competitive with each other; x° q polar (l^ max ) takes its maxima 
(minima) along q y = 0. This fact brings about a competition 
between the susceptibility and the inter-dimer interaction, and 
suppresses the instability toward the polar-CDW at low tem- 
perature. On the other hand, such competitive momentum de- 
pendences are weak for x q charge and V maK (q) near the non- 
polar CDW phase [see Figs [5J C ) an d ( e )L A similar reentrant 
feature and the competitions are suggested in Refs.l29landl33i 
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Finally, we mention the SC properties in the two-band 
model with t'/t < 1, that corresponds to k-(BEDT-TTF) 2 X 
for X=Cu(NCS) 2 and Cu[N(CN) 2 ]Cl. Phase diagram and gap 
functions calculated by RPA are presented in Fig. [6] In con- 
trast to the results in Fig. [2] where we take t'/t > 1, the SC 
phase with the d^^-type Bi 9 symmetry gap appears near 
the SDW phase. This change in the gap symmetry by chang- 
ing t'/t is consistent with the previous resultsJi^ The ex- 
tended s-type A% g gap is also realized near the CDW phase. 
As a consequence of the competition between the two SC 
phases with different symmetry gaps, a phase volume of the 
extended s-type SC is shrunk. We propose that the SC state 
induced by the polar charge fluctuation is favored in the case 
of t'/t > 1. 

In summary, we study the superconductivity in the low- 
dimensional organic conductor where electronic structure in- 
side of dimer is taken into account. We focus on roles of the 
polar charge fluctuation on the superconductivity. The two- 
types of the extended Hubbard models are analyzed by the 
RPA and FLEX methods. We find that the extended s-type SC 
state is realized near the CDW phase where the charge distri- 
butions inside of dimers are polarized. A monotonic increase 
of T c when the system approaches to the polar CDW phase 
is direct evidence of the SC state induced by the polar charge 



fluctuation. The obtained T c values are comparable to those in 
the k-(BEDT-TTF)2X systems. Competitive relationship be- 
tween this SC state and that induced by the spin fluctuation 
depends on the geometrical anisotropy in a triangle lattice; 
two SC states are cooperative (competitive) with each other 
in t'/t > 1 (t'/t < 1). The reentrant structures emerge in 
the boundary between the SC state and the polar-CDW state. 
The present study opens a new direction for research of the 
superconductivity and ferroelectricity. 
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